Landau Level Spectrum of ABA- and ABC-stacked Trilayer Graphene 
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We study the Landau level spectrum of ABA- and ABC-stacked trilayer graphene. We derive 
analytic low energy expressions for the spectrum, the validity of which is confirmed by comparison 
to a 7r-band tight-binding calculation of the density of states on the honeycomb lattice. We further 
study the effect of a perpendicular electric field on the spectrum, where a zero-energy plateau 
appears for ABC stacking order, due to the opening of a gap at the Dirac point, while the ABA- 
stacked trilayer graphene remains metallic. We discuss our results in the context of recent electronic 
transport experiments. Furthermore, we argue that the expressions obtained can be useful in the 
analysis of future measurements of cyclotron resonance of electrons and holes in trilayer graphene. 

PACS numbers: 81.05.ue, 71.70.Di, 73.43.Lp, 73.22.Pr 



I. INTRODUCTION 

Recent experimental realizations of graphene 
trilayersir— (TLG) have opened the possibility of 
exploring their intriguing electronic properties, which 
depend dramatically on the stacking sequence of the 
graphene layers X The low energy band structure for 
ABA-stacked TLG consists of one massless and two mas- 
sive subbands, similar to the spectrum of one single layer 
(SLG) and one bilayer graphenes (BLG), while ABC tri- 
layer presents approximately cubic bands.— Interestingly, 
when the TLG is subjected to a perpendicular electric 
field, a gap can be opened for ABC samplesj2i2i2rii 
similarly to bilayer graphene^ whereas ABA TLG 
remains metallic with a tunable band overlap J^ 

When a strong magnetic field is applied perpendicular 
to the TLG planes, the band structure is quantized into 
Landau levels (LLs) . The number of graphene layers as 
well as their relative orientation (stacking sequence) de- 
termine the features of the quantum Hall effect (QUE) 
in this material, where the Hall conductivity presents 
plateaus a t 14 ' 15 
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where N = 3 is the number of layers, n is the LL index, 
g = 4 is the LL degeneracy due to spin and valley degrees 
of freedom, — e is the electron charge and h is the Planck's 
constant. In particular, the plateau structure in a xy of 
TLG has been shown to be strongly dependent on the 
stacking sequenced 

In this paper we study the LL quantization of TLG. We 
obtain analytical expressions for the LL spectrum of TLG 
with ABA or ABC stacking order. The range of applica- 
bility of the analytical results is studied by a comparison 
to the density of states (DOS) obtained from a numeri- 
cal solution of the time-dependent Schrodinger equation 
within the framework of a tight-binding model on the 
honeycomb lattice* 1 ^— We further study the effect of a 
perpendicular electric field in the LL spectrum, finding 
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FIG. 1. Atomic structure of ABA- and ABC-stacked trilayer 
graphene. The intra-layer t and inter-layer 71 and 73 hopping 
amplitudes are schematically shown in the figure. 



that a zero-energy plateau develops in the Hall conductiv- 
ity only for ABC-stacked graphene, while ABA-stacked 
graphene remains ungapped. 

The paper is organized as follows. In Sec. |H] we ob- 
tain analytically the low energy LL spectrum of TLG. 
The analytic expressions of Sec. |H] are compared to the 
DOS numerically obtained from a full tight-binding cal- 
culation in the honeycomb lattice in Sec. IIIII Our main 
conclusions are summarized in Sec. IIVI 



II. 



ANALYTIC DERIVATION OF THE LANDAU 
LEVEL SPECTRUM 



In nature there are two known forms of stable stack- 
ing sequence in TLG, namely ABA (Bernal) and ABC 
(rhombohcdral) stacking. The difference between ABA 
and ABC stacking, schematically shown in Fig. [TJ is that 
the third layer is rotated with respect to the second layer 
by —120° (so that it will be exactly under the first layer) 
in ABA stacking, while it is rotated by +120° in ABC 
stackingiii^ In a basis with components of ipA ± , "0Si , 
Vu 2 ; ^b 2 j Vm 3 , ipB 3 , where ipA t (V's*) are the envelope 
functions associated with the probability amplitudes of 





tonian can be expressed in terms of a and a' as 



H 



( ° + 


A B a 











^ 


A B at 





7i 














7i 





A B a 





7i 








A B at 


























Asa 



V o 







71 



A 



, (3) 



Bfl 1 



J 



FIG. 2. (Color online) Low energy band structure of ABA- 
and ABC-stacked trilayer graphene around the K point. We 
have used the tight-binding parameters t — 3 eV and 71 = 
0.4 eV. The red dashed lines are a guide to the eye that mark, 
for the used parameters t and 71 , the position of the bottom 
(top) of the upper (lower) bands. The analytic expressions of 
these bands are given in Appendix A. 



the wave functions on the sublattice A (B) of the ith 
layer (i = 1,2,3), the effective low energy Hamiltonian 
of ABA-stacked TLG around the K point is£ 
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where p± = p x ± ip y , with p = (p x ,Py) the two- 
dimensional momentum operator, and vp = 3at/2 the 
Fermi velocity of the monolayer graphene, in terms of 
the in-plane nearest neighbor hopping t ra 3 eV and the 
carbon-carbon distance a w 1.42 A (from now on we use 
units such that h = 1 = c). For the moment, we only 
include the inter-layer hopping 71 w 0.4 eV in Eq. (J2J. 
The effective Hamiltonian for K' is obtained by exchang- 
ing p + and p- . The effect of far-distant hopping such as 
73 will be discussed in Appendix [Cj the Hamiltonian 
@ leads to a combination of two linear SLG-like bands 
[black lines in Fig. EJa)] and four massive BLG-like bands 
[red and green lines in Fig. Ufa)]. 

In the presence of an external perpendicular magnetic 
fieldj 2 ^ the canonical momentum p must be replaced by 
the gauge-invariant kinetic momentum p — > II = p + 
eA(r) where A(r) is the vector potential, and which obey 
the commutation relation [ILcIIy] = —i/l B , where Ib = 
\j\feB is the magnetic length. Therefore, this allows 
to introduce the ladder operators a = (?s/v / 2)n_ and 
a) = (/s/v / 2)n+, where W± = H x ±ili y) and which obey 
the commutation relation [0,0^] = 1. As in the usual 
one-dimensional harmonic oscillator, 



a \n) = \fn \n — 1) , a T \n) = y/n + 1 \n + 1) , 

where \n) is an eigenstate of the usual number operator 
a^a\n) = n\n), with n > an integer. Then, the Hamil- 



where Ab is the magnetic energy defined by As = 
V%vp/Ib- Therefore the six-components eigenstates 
of H can be reconstructed as ip = [cA 1 i Pn-i,kt 

i't 

CB 1 fn,k, CA 2 fn,k, CB 2 fn+l,k, CA 3 fn-l,k, CB 3 fn,k\ , 

where CAi(cBi) are amplitudes. If we choose the Lan- 
dau gauge A(r) = (0, Bx), then the wave function of the 
nth LL ip n ,k(x, y) is given by2i 

^*>ti-r{w^y'<*°-*^M> (4) 

where z — (x — kl B )/ls, H n (z) is the Hermitc polyno- 
mial, and ip n ,k = for n < 0. Then, the Hamiltonian 
matrix in the basis of ip is 
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with C\ — y/n and C 2 = \/n + 1. Eq. ([5]) has six eigen- 
values, which can be easily calculated: 



E n , s = ±^=[2 7l 2 + (2n+l)A 



+ Sv /4 7l 4 + 4 (2n + 1) 7l 2 A 2 B + A|] 1/2 , (6) 

E nfi = ±A B Vn~, (7) 

with s = ±1 and n > 0. The eigenstates correspond- 
ing to above LLs are given in Appendix B. Notice that 
Eq. ^ coincides (apart from a numerical factor v2 in 
front of 71) with the LL spectrum of a bilayer graphene, 22 
whereas the Eq. ([7]) corresponds to the LL spectrum of 
a single layer graphene. This is expected since the low 
energy band structure of ABA TLG consists of two mass- 
less SLG-like bands and four massive BLG-like bands, as 
it has been discussed above. In Fig. \5l&) we show the LL 
spectrum Eq. @-0 for ABA TLG obtained for the first 
50 LLs of each band (we only show the states with posi- 
tive energy). As in the zero magnetic field case, there are 
two sets of BLG-like LLs which disperse roughly linearly 
with B (the LLs plotted in red and green color) , whereas 
the linearly in k dispersing SLG-like band leads to a set 
of \AB-like LLs (plotted in black) [see Fig. (3{b) for a 
zoom of the low energy and low magnetic field region of 
Fig. [3] (a)]. Furthermore, a set of LL crossings occur due 
to the massless and massive characters of the subbands, 



as it has been observed experimentally^ Notice that the 
Landau levels in the low energy part of the spectrum have 
only E n _ character [see Fig. [3ja) and (b)], unless the 
magnetic field is very strong. For example, the third low 
energy Landau level belongs to the set of LLs -B n ,o when 
B > 45 T. On the other hand, the E n< + LLs only appear 

\E\ > |£o, + 



at an energy 
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In the limit 



nA 2 B <C 7 2 , the BLG-like bands Eq. ([6]) can be simplified 
to 
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which is similar to the commonly used expression for the 
low energy spectrum of BLG in a weak magnetic fields 
Whereas some of the results for the LL spectrum of 
ABA trilaycr graphene has been discussed before^ much 
less effort has been put on understanding the ABC TLG. 
However, recent experiments have shown the stability of 
TLG stacked with rhombohcdral order, and the possi- 
bility of opening a gap by applying a transverse electric 
field to the sample ^^ what has activated the interest 
on TLG with this stacking sequence. The Hamiltonian 
for ABC-stacked TLG around the K point is 
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The eigenvalues of Eq. ((5J) leads, as shown in Fig. [2] (b), 
to a low energy band structure that consists of a set of six 
cubic bands, two of them touching each other at the K 
point, and the other four crossing at an energy E = ±71 
above (below) the K point. In the following we will ob- 
tain the LL spectrum for this case. In a similar man- 
ner as for the ABA case, the six-components eigenstates 
of the Hamiltonian for ABC-stacked TLG can be recon- 
structed as if) = [cAiPn-l.k, CB 1 <Pn,k, CA 2 <£n,k, CB 2 <p n +l,k, 



CA 3 <£n+l,fc) CB i i Pn+2,k\ 

this case is (n > 0) 



and the Hamiltonian matrix in 
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(10) 
with C\ = y/n, C 2 = y/n+1 and C 3 = y/n + 2. The 
eigenvalues of Eq. ([TU]) are the solutions of the equation 

(11) 



E n + bE n + cE n + d = 0, 



where 

6=-2 7l 2 -3(l + 7i)A|, 
c = 1 l + 2(l + n) 1 2 1 A 2 B + 
d= -n(n+ l)(n + 2)A 



(2 + 6n + 3n 2 ) A%, (12) 
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which leads to a LL spectrum for ABC-stacked TLG 
given by2^ 
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E n3 =±i 2y/Qcos 



where 
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In Eq. (jlOp . the Landau level index n is required to 
be nonnegative. However, notice that Eq. (|10j) admits 
also eigenstates with real eigenvalues that contain com- 
ponents with n = — 1. The corresponding eigencnergies 
can be obtained by setting C\ = — 1, G-x = and C3 = 1 
in Eq. (|T0l) . This leads to three twofold eigenvalues that 
complement Eq. ([TB")) 

E-i,i=0, 



S_ li3 =±- v / 7 f + A|, 

where we label the contributions from the last two bands 
as -E-1,3, because they have a similar field dependence as 
the £„, 3 LLs [see Fig. G^d)]. 

In the low magnetic field limit, the Landau level spec- 
trum for ABC-stacked TLG can be approximated byi& 
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y/n(n + l)(n + 2). (17) 



The positive energy part of the LL spectrum obtained 
from Eq. p^|) is represented in Fig. Etc). One can dis- 
tinguish one set of LLs starting from zero energy, which 
correspond to the low energy band that touches the Dirac 
point, plus two set of LLs at an energy ~ 71 and which 
are related to the bands that cross at 71 [see Fig. |2fb)]. 
Whereas the low energy set of LLs can be understood 
from a standard quantization of a low energy cubic band, 
the LLs that appears at E n ~ 71 deserve some discus- 
sion [see Fig. GUd) for a zoom of the low field region of 
these states] . Most saliently, the hybridization of the up- 
per bands leads to two different sets of LLs. One set of 
LLs [plotted in green color in Fig. |3fc)-(d)], associated 
to the inner branches of the hybridized bands [denoted 
by the green lines in Fig. [UJb)], disperses with an energy 
E n > 71 and it is quite similar to that of a SLG. The 
other set of LLs, associated to the outer branches of the 
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FIG. 3. (Color online) Three band structures in the Landau level spectrum of the ABA- and ABC-stacked trilayer graphene. 
We have used Eq. (©-(0 f° r ABA stacking, and Eq. (|13p for ABC stacking. Only the first 50 Landau levels in each band are 
presented. 



hybridized bands [denoted by red lines in Fig. HJb)], has 
an energy that first decrease with B until it reaches a 
minimum value, and then grows in energy as B increases 
[see the lower set of LLs of Fig. [3Jd), which are colored 
in red] . This behavior is due to the cusp of this branch 
at £ = 71, and resembles the saddle point of the bilayer 
graphene bands in the presence of a transverse electric 
field. The effect of the perpendicular electric field in BLG 
is to open a gap in the spectrum, leading to Mexican hat 
like bands ) 12 i 25 ~ — with the corresponding anomalous LL 
quantization of the bandi 22 ' 31 ' 32 Therefore, the LLs as- 
sociated to the quantization of the lower branches of the 
hybridized bands in ABC TLG can be obtained, in a first 
approximation, by using the semiclassical approximation 
used in Ref. |3f| for a biased bilayer graphene. The de- 
generacy of zero-order Landau level in ABC TLG is three 
times larger than SLG. This result remains correct also 
for the case of inhomogenous magnetic field as follows 
from the index theorem^ 



III. DENSITY OF STATES FROM A FULL 
tt-BAND TIGHT-BINDING MODEL 

In order to check the range of validity of the ana- 
lytic expressions obtained in Sec. [Til in this section 
we compare the LLs obtained from the equations ([5])- 
([7]) and (|I3I) for the low energy spectrum of ABA- 
and ABC-stacked TLG, respectively, to the density of 
states (DOS) obtained numerically by solving the time- 
dependent Schrodinger equation (TDSE) on a honey- 
comb lattice in the framework of a 7r-band tight-binding 
modclJ£~— The effect of an external magnetic field is 
considered by means of a Peierls substitution 
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(18) 



where t mn is the hopping amplitude between sites m and 
n of the honeycomb lattice, and J A • d\ is the line in- 
tegral of the vector potential. A numerical study of the 
magneto-electronic properties of ABC TLG has been also 
reported in Ref. |3J. In Fig. [4] we compare our analytic 
results of Eqs. @-Q and §F^ with the numerical TDSE 
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FIG. 4. (Color online) Comparation of the Landau level spectrum obtained from the analytic expressions derived in the text 
(color lines) and the numerical simulation (black lines) of ABA- and ABC-stacked trilayer graphene. The sample used in the 
numerical simulations contains 3200 x 3200 atomic sites in each layer, and we use the periodic boundary conditions in the plane 
(XY) of graphene layers. 



results for the DOS, for two different values of magnetic 
field. We find a very good agreement between analytic 
and tight-binding results up to an energy of ~ 0.5 eV. 
Notice that when a LL crossing occurs, for example of a 
SLG-like LL crossing with a BLG-like LL in ABA TLG, 
this leads to an increase of the peak in the DOS. This is, 
e. g., the reason for the enhanced peaks at E ps 0.5 eV 
and E sa 0.7 eV in Fig. HJb), as it can be deduced by 
following the LL spectrum of Fig. [21a) at B — 50 T. 
Far from the neutral point, at an energy E > 0.5 eV the 
analytic results arc shifted to the right of the spectrum, 
as compared with the numerical TDSE results (see c. g. 
the peaks corresponding to £"„,- for ABA- and £"„,i for 
ABC-stacked TLG, represented by the red vertical lines 
in Fig. 0]). This is due to the fact that the dispersion 
relation for SLG is not linear anymore, so that higher or- 
der terms should be included for a precise reproduction 
of the position of the LLs. 

It is interesting also to check the range of validity of 
the most commonly used approximated expressions for 
the LL spectrum of TLG [Eq. © for ABA and Eq. $TQ 
for ABC]. Contrary to single layer graphene, for which 
the LL spectrum behaves as yBn up to rather high en- 
ergies (in Ref. |35| it was reported a deviation of only 
~ 40 meV at an energy of 1.25 eV), the B\Jn{n + 1) be- 
havior of the BLG-like LLs of ABA TLG as well as the 
B 3 / 2 y/n(n +l)(n + 2) behavior of ABC TLG are valid 
only in a rather reduced range of energies in the spec- 
trum. In fact, we see in Fig. [5] that, for the moderate 
value of magnetic field used for this plot (B = 20 T) the 



approximations Eqs. © and (fT?} fail to capture accu- 
rately even the second LL of the spectrum. The devia- 
tion is especially important for ABC trilayer graphene, 
as seen in Fig. [6l where one can see that there are devia- 
tions of hundreds of meV between the two results already 
for low LLs at some intermediate values of magnetic field 
~ 15 — 20 T. This is somehow expected since recent cy- 
clotron resonance experiment o 36 ! 37 on bilayer graphene 
required the use of the equivalent expression for BLG of 
Eq. ([5]), that we have obtained for the BLG-like bands of 
ABA TLG. Indeed, a good fitting (apart from some pos- 
sible many-body correction o 38 ' 39 ) of the magneto-optical 
experiments on BLG was achieved by using an expression 
similar to Eq. ©, with the only tight-binding parame- 
ters 70 = t and 71. Therefore, we expect that the ana- 
lytic expressions Eqs. ©-([7]) and (fT3| that we have ob- 
tained can be useful when analyzing future cyclotron res- 
onance experiments of ABA- and ABC-stacked trilayer 
graphene. 

Furthermore, motivated by recent transport measure- 
ments on TLG, which have revealed the strongly stack- 
ing dependent quantum Hall effect in this material^— we 
have calculated the Hall conductivity for the two stacking 
sequences of TLG, considering also the effect of a trans- 
verse electric field in the spectrum. Here the Hall con- 
ductivity <j xy is calculated by using the Kubo formula^ 
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where the charge density n s = L p (E) dE is obtained by 
integration of the DOS p(E) calculated from the TDSE 
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FIG. 5. (Color online) Landau level spectrum of (a) ABA- 
and (b) ABC-stacked trilayer graphene at B = 20 T obtained 
from the numerical solution of the TDSE using a 7r-band tight- 
binding model (black lines). In (a) the TDSE results is com- 
pared to the results from the analytic expression E n .- (red 
vertical lines) and E n ,o (green vertical lines) from Eq. ©-0! 
and to the approximation Eq. (0. In (b), the TDSE DOS 
are compared to the analytic result for E n ,i from Eq. 1)130 . 
and to the approximation Eq. (117)) (blue line). 



and 7r-band tight-binding method, and Aa xy is a correc- 
tion due to scattering of electrons with impurities^ and 
which is zero in the clean limit considered here. In Fig. 
we show the Hall conductivity of ABA- and ABC-stacked 
TLG with or without an external electric field. In the 
absence of any bias, the Hall conductivity for the two 
cases is similar, with plateaus at v = ±6, ±10, ±14, .... 
However, the structure of a xy is different when we con- 
sider the effect of a transverse electric field, which is ac- 
counted for here by adding a different (nonzero) on-site 
potential on the top and the bottom layers, namely, Ai/2 
on the top layer and — Ai/2 on the bottom layer. The 
main difference between ABA- and ABC-stacked TLG 
in the presence of a transverse bias is that it leads to 
a gap opening in the case of ABC-stacking, while the 
ABA-stackcd TLG remains gapless, as it has been ob- 
served experimentally^ In fact, the opening of the gap 
and the corresponding insulating state leads to the ap- 
pearance of a zero energy plateau in the Hall conductiv- 
ity in ABC TLG, plateau which is absent in ABA TLG, 
as shown in Fig. [7] for different values of Ai. On the 
other hand, the position of the plateaus depends very 
much on the value of the induced difference potential 
Ai. For a small bias leading to Ai = 0.15 eV, we find 
plateaus for ABA TLG at v = ±2, ±4, ±6, ±10, ±14,..., 
whereas a higher value, Ai = 0.3 eV leads to plateaus 
at v = ±2, ±6, ±8, ±12, ±14,.... On the other hand, 
whereas Ai = 0.15 eV leads to plateaus for ABC at all 
even values of v (including v = 0), some of the plateaus 
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FIG. 6. (Color online) Comparation of the analytic results 
for the first ten Landau levels obtained from Eqs. ©-© an d 
(|13)l for ABA- and ABC-stacked graphene, respectively (black 
solid lines) and the approximations Eqs. (jSJ) and (|17p (red 
dashed lines). The inset panels are the difference between 
Eq. ()6)l- (|13)) and the commonly used approximations Eqs. 
l[S)l-(|17p. Notice the different range of magnetic fields used in 
the inset with respect to the main figures. 



are missing for a higher value of bias, Ai = 0.3 eV, for 
which we find plateaus at v = 0, ±2, ±4, ±6, ±12, ±16, .... 
In fact, a more deep understanding of the Hall conduc- 
tivity of TLG would require further analysis, which is 
beyond the scope of this work. Furthermore, we empha- 
size that even experimentally, there is no consensus so 
far about the structure of the quantum Hall plateaus in 
trilayer graphene, having been found different structures 
for almost every transport measurement £2r— 



IV. CONCLUSIONS 

In conclusion, we have derived analytic expressions for 
the Landau level spectrum of trilayer graphene. The two 
stable stacking sequence, ABA (Bcrnal) and ABC (rhom- 
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FIG. 7. (Color online) Hall conductivity of ABA- and ABC- 
stacked trilayer graphene with different values of Ai induced 
by a transverse electric field. 



equation of a 7r-band tight-binding model on the honey- 
comb lattice. We find very good agreement between the 
numerical solution and the analytic approximation for 
the spectrum up to an energy of ~ 500 meV. However, 
we show that the most commonly used approximations 
for the spectrum of TLG, for which th e BLG-like LLs 
of ABA TLG disperse as By/n(n + 1) and the LLs for 
TLG disperse as B 3 l 2 ^/n(n + l)(n + 2), fail to capture 
even the lower LLs already for moderate magnetic fields 
of ~ 20 T. Therefore, we believe that our results maybe 
useful for the analysis of future magneto-optical measure- 
ments, which has been successfully applied to study the 
LL spectrum of SLGii^ and BLG&2L 

Finally, we have calculated the Hall conductivity of 
TLG by means of the Kubo formula. The inclusion of a 
transverse electric field leads to a gap opening in ABC 
TLG, whereas ABA TLG remains metallic. This effect 
is seen by the appearance of a zero energy plateau only 
for ABC stacking, in agreement with recent transport 
experiments ££~— 
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Appendix A: Band structure of ABA and ABC 
trilayer graphene in the absence of magnetic field 



bohcdral) have been considered. The LL spectrum for 
ABA TLG is composed by a set of bilayer graphene-like 
LLs, which disperse at low energies as B, and a set of sin- 
gle layer graphene-like LLs, which disperse as y/~B. The 
different character of the bands lead to a series of LL 
crossings, which has been observed experimentally^ On 
the other hand, the six cubic bands of ABC TLG leads 
to a rather peculiar LL quantization of the spectrum. 
Whereas the bands that touch the Dirac point lead to a 
set of B 3 ! 2 LLs, the hybridization between the two bands 
that cross each other at E = 71 leads to one set of mass- 
less like LLs (with energy E > 71 ) , and a set of LLs which 
present a minimum and then grows with B, associated to 
the lower branch of the hybridized bands. The presence 
of the minimum on this set of LLs is associated to the 
presence of a cusp in this branch of the spectrum, in a 
similar manner as the Mexican hat like dispersion of a 
biased bilayer graphene. 

The range of validity of our analytical results is checked 
by comparing the LL spectrum obtained in the contin- 
uum approximation to the density of states obtained from 
the numerical solution of the time dependent Schrodingcr 



In the absence of a magnetic field, the Hamiltonian of 
ABA-stacked TLG around the K point is given in Eq. 
@, with eigenenergies given by 



E s = ±[ 7l 2 + v 2 F k 2 + aJ-rt + 2j 2 v 2 F k 2 } 1/2 , s = ±1 



Eq = ±VFk. 



(Al) 



Similarly, for ABC-stacked TLG, the Hamiltonian Eq. 
(O leads to the eigenvalue problem 

E 6 - (2 7l 2 + 3v% k 2 ) E A + 

( 7l + 2 1 2 v 2 F k 2 + 3v F k 4 ) E 2 - v 6 F k 6 = 0, (A2) 

the solutions of which take the form of Eq. (fTB"]) with the 
new quantities b = —2-/ 2 — 3v F k 2 , c = yf + 2~/ 2 v F k 2 + 
3v F k 4 and d = —v F k 6 . In fact, Eq. (|A2[) can be decom- 
posed into the two equations: 



E 3 + v F kE 2 



7i 



v 2 F k 2 ) E - v 3 F k 3 



0, (A3) 



E s 



v F kE 2 - ( 7l 2 + v 2 F k 2 ) E + v 3 F k 3 = 



(A4) 



the solutions of which arc 



where s = ±1 correspond to the solutions of Eq. (|A3 
and (|A4[) respectively, in terms of the new parameters 



6 = cos" 1 



sR 



(A6) 



E a . = 2^JQco& 



4-27r\ vpk 



I— f9 + 4n\ v F k 
Ep ja =2VQcos I — - — I -a-g-. 

„ /-r f6\ V F k 

E-, >a = 2 v /Qcosl - I -s-g- , 



(A5) 



i? = 



8v F k 3 vpk"f\ 
37? + 4v 2 F k 2 



(A7) 
(A8) 



Appendix B: Wave functions of ABA trilayer 
graphene 

From the matrix Hamiltonian Eq. ([3]) one can calculate 
the eigenstates of the ABA TLG. They are given by 



tpn, s {x,y) = 



± 



nAi 



\ATB„ iS Ai; 



En. 



-1 



and 



1 



{l+n)A%(n A B 2 ~El s ) 
7l 2 E 2 ~s 



y/nA 

(l+n)A 2 B (n A B 2 -El s ) 



</' 



± 



71 -c/„ 



± 



± 
A B 2 -B, 2 



-nAi 



lAt 



7i 7i ^ 

/1+iiAb _ nyi+nA B 3 
71 E 2 s 



71 



71 
Vn,k( X >y) 



71 
¥>n,k( X 'V) 



tp n -i t k(x,y) 



71,0 



± 



-E„ , 



^„+i.fc(a;, y) 



<^«+2,fe(a;,y) 



(Bl) 



Tfn-l,k{ X >y) 

-Vn,k{x,y) 




±¥> n+ i,k(a;,y) 
Vn+2,k(x,y) 



(B2) 



Notice that the states with the eigenvalues E n> o are 
the surface states which arc located only on the top and 
bottom layers, and these surface states in each layer have 
the same expressions as the single-layer graphene. 



Appendix C: Effect of 73 in the DOS 

In this appendix we study the effect of considering, 
besides t and 71 , the inter- layer hopping amplitude 73 in 
the spectrum (see Fig. [1]). In Fig. [51 we compare the 
Landau level spectrum and Hall conductivity of ABA- 



r 



and ABC-stacked TLG with and without 73. Here we 
use 73 = 0.3 cV as it is in the nature graphite i 43 i 44 For 
the considered magnetic field, the effect of 73 in the spec- 
trum is negligible, as seen in Fig. [5] Therefore, trigonal 
warping has very small effect to the low energy spectrum 
of the Landau levels in the presence of high magnetic 
field. In fact, this is also the case in bilayer graphene, 
where the LL spectrum can be adequately described by 
neglecting 73 over the field range where l^ 1 > §073771 
(where m ~ 0.054m e is the effective mass in the bulk 
graphite)^ In our calculations, the DOS and Hall con- 
ductivity are almost the same, as it is shown in Fig. [5] 
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FIG. 8. (Color online) Comparison of the Landau level spectrum and Hall conductivities of ABA- and ABC-stacked trilayer 
graphene with (red lines) or without (black lines) considering the interlayer hopping parameter 73. 



2 W. Bao, L. Jing, Y. Lee, J. V. Jr., P. Kratz, D. Tran, 
B. Standley, M. Aykol, S. B. Cronin, D. Smirnov, 
M. Koshino, E. McCann, M. Bockrath, and C. Lau, 
(2011). larXiv:1103.6088l 

3 C. H. Lui, Z. Li, K. F. Mak, E . Cappelluti, and T. F. 
Heinz, (2011). larXlv:1105.4658l 

4 A. Kumar, W. Escomer, J. M. Poumirol, C. Faugeras, D. P. 
Arovas, M. M. Fogler, F. Guinea, S. Roche, M. Goiran, 
and B. Raquet, (2011'). larXiv:1104.1020"1 

5 L. Zhang, Y. Zhang, J. Camacho, M. Khodas, and I. A. 
Zaliznyak, f20in, larXiv:1103.6023l 

6 S. H. Jhang, M. F. Craciun, S. Schmidmeier, S. Toku- 
mitsu, S. Russo, M. Yamamoto, Y. Skourski, J. Wos- 
nitza, S. Tarucha, J. Eroms, and C. Strunk, (2011), 
larXiv:1106.4995l 

7 F. Guinea. A. H. Castro Neto, and N. M. R. Peres, Phys. 
Rev. B 73, 245426 (2006). 

8 M. Koshino and E. McCann, Phys. Rev. B 80, 165409 
(2009). 

9 A. A. Avetisyan, B. Partoens, and F. M. Peeters, 



Phys. Rev. B 79, 035 421 (2009) 
A" A" AVetisyan, B. Partoens, 



and F. M. Peeters, 



Phys. Rev. B 81, 115432 (2010) 



F. Zhang, B. Sahu, H. Min, and A. H. MacDonald 



Phys. Rev. B 82, 035409 (2010) 

E. V. Castro, K. S. Novoselov, S. V. Morozov, N. M. R. 

Peres, J. M. B. L. dos Santos, J. Nilsson, F. Guinea, A. K. 

Geim, and A. H. C. Neto, Phys. Rev. Lett. 99, 216802 

(2007). 

M. F. Craciun, S. Russo, M. Yamamoto, J. B. Oostinga, 

A. F. Morpurgo, and S. Thrucha, Nat. Nanotechnol. 4, 

383 (2009). 



14 E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96, 086805 
(2006). 



H. 



Min 



and 



A. 



H. 



MacDonald, 



16 

17 S 



Phys. Rev. B 77, 155416 (2008) 



18 



A. Hams and H. De Raedt, |Phys. Rev. E 62, 4365 (2000) 
Yuan, H. De Raedt, and M. I. Katsnelson, 
|Phys. Rev. B 82, 115448 (2010)[ 

S. Yuan, H. De Raedt, and M. I. Katsnelson, 
|Phys Rev. B 82, 2 35409 (2010) 

19 M. Koshino, Phys. Rev. B 81, 125304 (2010). 

20 For a recent review on the electronic properties of 
graphene in a magnetic field, see. M. O. Goerbig, (2010), 
larXiv: 1004.33961 

21 M. Koshino and T. Ando, Phys. Rev. B 77, 115313 (2008). 

22 J. M. Pereira, F. M. Peeters, and P. Vasilopoulos, Phys. 
Rev. B 76, 115419 (2007). 

23 M Koshino and E. McCann, 

Phys. Rev. B 83, 165443 (2011) 1 



24 See e.g. , G. Birkhoff, and S. M. Lane, A Survey of Modern 
Algebra, 5th ed. (New York: Macmillan, 1996). 

25 E. McCann, |Phys. Rev. B 74, 161403 (2006) | 

26 H. Min, B. Sahu, S. K. Banerjee , and A. H. MacDonald, 
|Phys. Rev. B 75, 155115 (2007)1 

27 Y. Zhang, T.-T. Tang, C. Girit, Z. Hao, M. C. Martin, 
A. Zettl, M. F. Cromm ie, Y. R. Shen, and F. Wang, 
|Nature 459, 820 (2009)] 

28 K. F. Mak, C. H. Lui, J. Shan and T. F. Heinz, 
|Phys. Rev. Lett. 102, 256405 (2009)] 

29 A. B. Kuzmenko, E. van Heumen, D. van der Marel, 
P. Lerch, P. Blake, K. S. Novoselov, and A. K. Geim, 
|Phys. Rev. B 79, 115441 (2009)| 

30 D. Wang and G. Jin, Euro. Phys. Lett. 92, 57008 (2010). 



10 



31 L. M. Zhang, M. M. Fogler, and D. P. Arovas, (2010), 
larXiv: 10 08. 1418 



M. Koshino, (2011). larXiv7l l05. 5919 



M. 



1. 



Katsnelson and 



M. 



Prokhorova, 



|Phys. Rev. B 77, 205424 (2008) 

34 C. Ho, Y. Ho, Y. Chiu, Y. C hen, and M. Lin, 
| Annals of Phy sics 326 , 721 (2011)| 

35 P. Plochocka, C. Faugeras, M. Orlita, M. L. Sad- 
owski, G. Martinez, M. Potemski, M. O. Goer- 
big, J.-N. Fuchs, C. Berger, an d W. A. de Heer, 
|Phys. Rev. Lett. 100, 087401 (2008)] 

36 E. A. Henriksen, Z. Jiang, L.-C. Tung, M. E. Schwartz, 
M. Takita, Y.-J. Wang, P. Kim, and H. L. Stormer, Phys. 
Rev. Lett. 100, 087403 (2008). 

37 M. Orlita, C. Faugeras, J. Borysiuk, J. M. Baranowski, 
W. Strupihski, M. Sprinkle, C. Berger, W. A. de Heer, 



D. M. Basko, G. Martinez, 
B 83, 125302 (2011). 
38 R. Roldan, J.-N. Fuchs 



and M. Potemski, Phys. Rev. 
and M. O. Goerbig, 



40 



Phys. Rev. B 82, 205418 (2010) 



K. Shizuya, (2011), arXiv:1103.5696 

R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 

M. L. Sadowski, G. Martinez, M. Potemski, C. Berger, and 

W. A. de Heer, Phys. Rev. Lett. 97, 266405 (2006). 

Z. Jiang, E. A. Henriksen, L. C. Tung, Y.-J. Wang, M. E. 

Schwartz, M. Y. Han, P. Kim, and H. L. Stormer, Phys. 

Rev. Lett. 98, 197403 (2007). 

B. Partoens and F. M. Peeters, 



Phys. Rev. B 74, 075404 (2006) 

A" EL Castro-Neto, F. Guinea, N. M. R. Peres, 

K. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009). 



